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ABSTRACT 


Steam leakage from an aircraft carrier catapult is sometimes ingested by the 
aircraft’s engines upon launch which may induce compressor stall. Investigation of the 
phenomenon known as a “pop stall” is of particular importance as the Navy prepares to 
field the F35C, the aircraft carrier variant of the joint strike fighter. The single engine 
design of the F-35C makes this aircraft particularly susceptible to steam-induced stall 
during catapult launch. The present project examined compressor stall and included 
steady-state as well as transient measurements in the inlet of a transonic compressor prior 
to and during a steam-induced stall. Hotwire measurements of the inlet flow field were 
taken to determine an inlet turbulence intensity of 2-3% during both subsonic as well as 
transonic compressor operation. A 95% speed line was established from data taken from 
open throttle to near stall. Hot-film and Kulite pressure data taken near stall showed the 
existence of a stall precursor which appeared near half rotor speed. Steam was injected 
into the inlet; however the initial method added mass to the system and did not induce a 
stall. A decrease in the amplitude of the pressure trace was observed however. A stall 
was induced by steam ingestion ahead of the existing inlet throttle, with upstream 
transient measurements taken using both hot-film and Kulite pressure transducers. 
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I. INTRODUCTION 


The catapult-launch systems on today’s aircraft carriers are known to release 
steam during the launch cycle, especially as the seals degrade over time. When this 
steam is ingested into the intakes of aircraft, it can cause a phenomenon known as “pop 
stall” causing a loss of power and possible total engine stall. Experiments with an F-18 
conducted at Lakehurst Naval Engineering Station (Figure 1) demonstrated the 
susceptibility of current aircraft to steam induced stall. 



Figure 1. F-18 “Pop-Stall” during Fakehurst steam ingestion experiment 

This phenomenon is of particular concern to the Navy as it prepares to transition 
to the single engine F-35C, the carrier variant of the Joint Strike Fighter (JSF). The 
single engine design of this aircraft increases the probability of a pop stall resulting in 
catastrophic loss of the aircraft. 

The focus of the work conducted at the Turbopropulsion Faboratory (TPF) at the 
Naval Postgraduate School (NPS) is concentrated on the “pop-stall” problem. 
Investigations conducted with the Transonic Compressor Rig (TCR) are intended to 
improve the understanding of the nature of steam induced stall. The present transonic 
compressor fan stage was designed specifically for the Naval Postgraduate School to be 
used in their TCR by Sanger (1996) at the NASA Glenn Research Center (Ref 1 and 2). 
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Extensive amounts of data were eollected by Gannon, Hobson and Shreeve in 
order to map the performance characteristics of the compressor in both the fan-stage as 
well as rotor-only configurations. The data were used to investigate the conditions of the 
TCR prior to and during stall as well as establish performance characteristics at 70%, 
80%, 90% and 100% design speed. (Ref 3-5) 

Unsteady pressure measurements at 60%, 70%, and 80% design speed were 
reestablished by Rodgers in 2003 (Ref 6). Villescas conducted inlet and exit surveys at 
70%, 80%, 90% and 100% design speed with a three-hole probe and determined 
spanwise distributions of the rotor diffusion factor at choke, peak efficiency and stall 
(Ref 7). Brunner repeated Villescas’ inlet flow-field surveys with a five-hole probe and 
determined pitch angle and Mach number distributions in the inlet to the rotor (Ref 8). 

The present experiment furthered the understanding of the performance 
characteristics of the transonic compressor rotor prior to and during steam-induced stall. 
Turbulence intensity measurements were taken upstream of the rotor in order to 
characterize the inlet flow conditions, as well as to document the inlet turbulence for 
future computational fluid dynamics (CFD) simulations. Performance measurements 
were taken at 95 percent speed from open throttle to near stall, and the data were added to 
the existing compressor map. Data taken upstream of the rotor at 90% speed by both a 
hot-film probe as well as Kulite pressure transducers in the case wall were analyzed to 
determine stall precursor events. Finally, transient data were collected from both the hot- 
film as well as Kulite pressure transducers during steam induced stall at 70% speed. 
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II. EXPERIMENTAL FACILITY AND EQUIPMENT 


The present study was condueted at the Turbopropulsion Laboratory (TPL) within 
the Department of Meohanieal and Astronautieal Engineering a the Naval Postgraduate 
Sehool (NPS). 

A. TRANSONIC COMPRESSOR 

The Transonic Compressor Rig (TCR), as shown in Figure 2, was driven by two 
opposed-rotor air turbine stages, supplied by an Allis-Chalmers axial compressor. The 
Allis-Chalmers compressor supplied three atmospheres of air pressure at a mass flow rate 
of 5 kg/s. The configuration tested was detailed extensively by O’Brien (Ref. 9) and 
Papamarkos (Ref 10). 



Figure 2. Transonic compressor rig in test cell with inlet piping removed. 

The transonic compressor rotor, was designed as part of a complete fan stage by 
Sanger using a CFD code (Ref 1). By employing numerical machining, the stage was 
manufactured specifically for testing and evaluation within the TCR. 
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Figure 3. As tested rotor only eonfiguration. 


The rotor had 22 blades and was made from a high strength aluminum alloy 
(7075-T6). In the present experiment, the rotor was tested with a parabolie spinner, 
whieh replaeed the eonieal spinner used by O’Brien (Ref. 9). 

In most of the previous studies the entire stage was evaluated. In the present 
experiments, the stator was removed and only the rotor was present. This was to use the 
simplest eonfiguration during the initial steam ingestion experiments. Table 1 below was 
reprodueed from Sanger’s “Design Methodology” (Ref 1) and details the design 
parameters for the rotor. 
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Table 1. Sanger rotor design parameters (Ref 1) 
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B. STEAM INGESTION/COMPRESSOR SYSTEM 

The steam generator used in the present experiment, as seen in Figure 4, was a 
Sussman model SVS600 steam boiler. The SVS600 was capable of producing saturated 
steam up to a maximum working pressure of 1000 kPa or 1.4 kg/sec at lOOC. (Ref 11) 
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Figure 4. Sussman model SVS600 steam boiler. 
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Figure 5. Steam ingestion/compressor system 
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The layout of the compressor rig and steam ingestion system is shown in Figure 5. 
Steam produced by the steam generator was directed via a 7.62 cm(three inch) diameter 
pipe and vented just outside the intake plenum. Figure 6 shows the orientation of the 
steam pipe with respect to the intake plenum. A pressure transducer was installed in the 
steam pipe so that the transient response of the steam pressure could be monitored. Two 
fast-acting, remote-operated solenoid valves were used for venting the pipe as well as 
releasing steam into the intake plenum. 



Figure 6. Steam pipe and intake plenum orientation 


Air was drawn into the compressor from atmosphere through a hydraulically 
operated throttle valve shown in Figure 5. A five-meter long 46cm diameter pipe 
connected the settling chamber to the test compressor. The pipe contained a nozzle, 
which was used for flow rate measurements, as well as a transition duct to the 27.94cm 
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case wall. A hot-film probe and a thermoeouple probe were mounted in the case wall 
11.5 cm upstream on the leading edge of the rotor, as shown in Figure 5. 
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Ill INSTRUMENTATION 


A. PROBES 

1. Hot-Film 

Hot-wire anemometry measurements were taken in order to characterize the inlet- 
flow conditions. Measurements were taken using a 20 micron TSI model 1212-20 hot- 
fllm probe (Figure 8, Ref 12) with serial number 70536123. Data were acquired using an 
IFA 100 Intelligent Flow Analyzer, and TSI ThermalPro software version 4.53, running 
on a PC. 


36 mm {I.SOJ 


rC 


7.51 mm (.3.12)1 
- 1 . 


LJQ 


mm (.125) Dia. ~ mm (.16) Dia. 


Figure 7. TSI model 1212-20 hot-fllm probe (From Ref 12) 

2. Kulite Pressure Transducer 

A Kulite Pressure Transducer, model XCQ-080, was used to obtain time-resolved 
wall pressure data. The probe was a miniature, semiconductor, strain-gauge transducer 
that incorporated a fully active four-arm Wheatstone bridge, which was dielectrically 
isolated on a silicon-on-silicon diaphragm. The natural frequency of the probe was 300 
kHz. A diagram of the probe is given in Figure 8 and the factory speciflcations are given 
in Table 2. 
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Figure 8. Kulite XCQ-080 series transducer (Ref 6) 


Table 2. Factory specifications for XCQ-080 (Ref 6). 


Input 

Pressure Range 

25psi 

Over Pressure 

50 psi 

Burst 

75psi 

Rated Electrical Excitation 

10VDC/AC 

Maximum Electrical Excitation 

15VDC/AC 

Input Impedance 

800 Ohms 

Output 

Output Impedance 

1000 Ohms 

Full Scale Output 

100 mV 

Residual Unbalance 

+-3% FSO 

Non-Linearity and Hysterisis 

0.1 %FS BFSL 

Hysteresis 

0.1% 

Repeatability 

0.1% 

Resolution 

Infinite 

Natural Frequency 

300kHz 

Perpendicular Accel Sensitivity 

0.0003% FS/g 

Transverse Accel Sensitivity 

0.00004% FS/g 

Insulation Resistance 

100 Megohm 

Environmental 

Operating Temp Range 

-65 to 250 deg F 

Compensated Temp Range 

80 to 180 deg F 

Thermal Zero Shift 

+-1 % FS/100 F 

Thermal Sensitivity Shift 

+-1 % FS/100 F 
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B. 


INSTALLATION OF SENSORS 


1, Hot-Film Probe 

The hot-film probe was installed in the ease wall 11.5em upstream of the leading 
edge of the rotor, and extended 4cm into the flow as shown in Figure 5. As can be seen 
in Figure 7, the sensor needles that hold the hot-film were bent forward and in this 
experiment the probe was positioned so that the leads pointed into the flow. The probe 
holder used to secure the probe is shown in Figure 9. The relative position of the probe 
ahead of the rotor in the inlet casing can be seen in Figure 10. 



Figure 9. Flot-film probe holder installed in inlet casing 
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Figure 11. Relative positions of Kulite pressure transdueer and blades. (From Ref 6) 


Figure 10. Relative position of hot-film probe 


2, Kulite Pressure Transducers 

Three Kulite pressure transdueers (designated 2, 8, and 9 in Figures 11 and 12) 
were mounted flush with the ease wall. Transducer 2 was located closest to the rotor 
approximately 1 cm from the leading edge just off top dead center. Transducers 8 and 9 
were located 120 degrees from transducer 1 but approximately 2 cm ahead of the rotor. 
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Figure 12. Relative positions or Kulite pressure transdueers in ease wall (From Ref 6) 

C. DATA ACQUISITION 

1. Hot-Film Data Acquisition 

Figure 13 shows the layout of the hot-film data aequisition system. An IFA 100 
Intelligent Flow Analyzer was used for hot-film data aequisition. The probe was 
eonneeted to the IFA 100 via a double-shielded eoax eable. Instantaneous voltages from 
the anemometer were sent to a United Eleetronie Industries Power DAQ A/D Board 
Model PD2-MFS-4-1M/12 installed in the P.C. via a terminal board. Control of the IFA 
100 from the PC was via an RS-232-C Digital Control Line. 
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Figure 13. Hot-film data acquisition system 

2. Kulite Data Acquisition 

The Kulite data acquisition system is shown in Figure 14. The Kulite pressure 
transducer was connected to the Hewlett-Packard E1529A Remote Strain Conditioning 
Unit via an RJ-45 LAN cable. Also connected to the Hewlett-Packard E1529A was a two 
channel power supply used to meet the 5 volt bridge excitation requirement of the 
Kulites. The HP E1529A was connected to the data interface port of the HP E1422A via 
an RJ-45 cable. The HP E1529A was connected to the HP E1433A digitizer via a 37-pin 
connection. A tachometer signal was connected to the HP E1433A in order to provide a 
speed reference to the data. Also connected to the HP E1433A was the signal from the 
lEA 100. This allowed for easy time-based comparison of the hot-film data and the 
Kulite data. Both the lEA 100 signal as well as the tachometer signal came in via 
standard coax cables and a break-out box was used as an adapter to route these signals 
into the 27 pin connection of the HP El433A. The HP El433A and HP El422A were 
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addressed through the HP E8404A VXI Mainframe and interfaced to a PC. A full 
description of the Kulite data acquisition system was given by Rodgers (Ref 6). 
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Figure 14. Kulite data acquisition system 
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IV. EXPERIMENTAL PROCEDURE 


A. PROBE AND TRANSDUCER CALIBRATION 

Calibration of the hot-film probe was eondueted in situ. After inserting a shorting 
probe in the plaee of the hot film probe, the eable resistanee was measured using the IFA 
100. After replaeing the shorting probe with the hot-film probe, the probe resistanee was 
measured. Operating resistanee and bridge eompensation was set aeeording to 
manufaeturer speeifieations, whieh were probe speeifie. The low and high flow voltage 
was measured with the eompressor running at 30% andl00% speed. The eompressor 
speed was then returned to 30% and inereased inerementally as measurements were 
taken. With eaeh hot-film measurement, the static to stagnation pressure difference was 
entered into the ThermalPro software for a velocity calculation. After all of the 
calibration measurements were taken, the ThermalPro software applied a curve fit to the 
calibration data and King’s Law coefficients were calculated. A complete description of 
the calibration procedures can be found in Appendix A. 

Calibration of the Kulite was carried out while the compressor was running, by 
applying different reference pressures to the reference tube, (Figure 8) and averaging the 
voltage recorded by the data acquisition system. The temperature dependence of the 
Kulite was alleviated by calibrating while online. A detailed account of the calibration 
procedure is given by Rodgers (Ref 6). 

B, TYPICAL COMPRESSOR OPERATION 

Standard operating procedure for testing the TCR was to maintain a constant rotor 
speed while taking measurements at different throttle settings. By closing the throttle, the 
mass flow rate was reduced and the rig operating point could be determined by the 
procedures described by Gannon et al. (Ref 6). 

The first step was to open the throttle and adjust the rotor RPM to the desired 
speed. Table 3 is a list of typical rotor speed settings. Measurements were then taken 
with the rotor speed constant and the mass flow rate varied by actuating the hydraulic 
throttle (Figure 5). Mass flow rate, inlet and exit total temperatures and pressures were 
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measured to ealculate total-to-total pressure ratio and isentropie efficiency in order to 
determine the position on the compressor performance map. Measurement and 
calibration procedures were described in more detail by Gannon et al. (Ref. 3). Along 
each speed line, measurements were taken from full open throttle to near stall. 


Tables. TCR speed settings 


RPM 

Percent Speed 

27,085 

100 

25,730 

95 

24,376.5 

90 

21,668 

80 

18,959.5 

70 


C. HOT-FILM MEASUREMENTS AT 70, 90, AND 95 PERCENT SPEED 

Data at 70, 90, and 95 percent speed were obtained in order to characterize the 
turbulence intensity levels of the inlet flow field. Data were recorded at a sampling rate 
of 100,000 Hz from the probe. The duration of each measurement sample was 0.0410 
seconds, and contained 4096 values. Three data samples were collected at 70 percent 
speed at open throttle, peak efficiency, and near stall. At 90 and 95 percent speed, data 
were taken over the whole operating range from open throttle to near stall. 

At 90 percent speed, several measurements were also made with the hot-film 
probe for a longer sample time. Specifically, the sample duration was 5.2429 seconds 
and consisted of 524,288 data points. Full procedures for operation of the ThermalPro 
software are given by Brown (Ref 13). 

D, PERFORMANCE MEASUREMENTS AT 95% SPEED 

The procedure for creating the 95 percent speed line was similar to that mentioned 
above. The compressor was brought to 95 percent speed with the throttle completely 
open. Incrementally, over 16 points, the throttle was closed and mass flow rate, total to 
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total pressure ratio, and isentropic effieiency were measured using a proeedure described 
by Gannon et al. (Ref 3). The compressor was throttled to a point near the expected stall 
point but not actually stalled in order to prevent damage to the rotor. 


E. STEAM-INDUCED STALL RUNS 

Measurements were conducted in order to ascertain the point at which steam 
ingested into the compressor would induce stall. For this part of the experiment, the 
compressor was initially set at 70 percent speed and peak efficiency. The throttle was 
then closed incrementally as measurements were taken at each setting. For each throttle 
setting, an initial steady-state measurement was taken prior to steam ingestion. Again, 
mass flow rate, total-to-total pressure ratio, and isentropic efficiency were measured in 
order to determine the location of the setting on the compressor map. Hot-film data 
samples lasting 5.2429 seconds, that consisted of 524,288 data points, were also taken. 
Kulite data were also collected at a rate of 196,608 Hz at sample lengths that varied from 
3 to 8 seconds. 

For each throttle setting, and after the steady-state measurements were taken, the 
following procedure was followed for steam ingestion. The boiler isolation valve and 
steam vent solenoid valve were both opened to allow the steam pipe to warm up. After 
the pipe was fully heated, the vent valve was closed and the data trace from the pressure 
transducer and thermocouple was started. When the pressure in the pipe reached a target 
pressure, the isolation valve was closed. A three second countdown commenced. At two 
seconds prior to start, the Kulite data acquisition was initiated. At one second prior to 
start, the hot-film data acquisition commenced. At the end of the countdown, the main 
steam valve (fast-acting solenoid valve) was opened allowing the steam to dump into the 
inlet plenum of the compressor. The steam was (observation showed) entrained through 
the throttle into the inlet by the suction of the compressor. After several seconds, all data 
acquisition was stopped. The whole procedure was repeated at reduced throttle settings 
until stall occurred when the steam was ingested by the compressor. 

Post processing of the data was conducted with MATLAB (Ref. 14). Specific M- 
files that were used are given in Appendix C. 
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V. RESULTS AND DISCUSSION 


A, TURBULENCE INTENSITY MEASUREMENTS 

Turbulence intensity values were calculated from hot-film data taken at 70, 90, 
and 95 percent speeds. These values were then plotted versus mass flow rate as seen in 
Figure 15. 



Figure 15. Turbulence intensity vs. mass flow rate 

An average turbulence intensity of 2.26 percent with a variance of 0.287 was 
determined from these measurements. There does not appear to be a correlation between 
turbulence levels, and mass flow rate or speed. 
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B, 95% SPEED PERFORMANCE DATA 

Extensive testing of the eompressor was eonducted by Gannon et al. at 70, 80, 90, 
and 100 percent speeds (Ref 3). It was eonsidered desirable to obtain performance data at 
95 percent speed as it was comparable to cruise conditions for an actual engine. Figures 
16 and 17 show the compressor performance map with the addition of the 95 percent 
speed line. These data closely follow those previously taken at 90 and 100 percent speed 
both with respect to the trend of the total pressure ratio characteristic and the value of the 
measured peak efficiency of 88.6 percent at 7.5kg/sec mass flow rate. The peak total-to- 
total pressure rise at the 95 percent speed was measured to be 1.56 at a mass flow rate of 
6.5 kg/sec 


22 




Figure 16. Pressure ratio versus mass flow rate with 95% speed line. 
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Figure 17. Isentropie efficieney versus mass flow with 95% speed line. 
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C. MEASUREMENTS OF STALL PRECURSOR WITH A HOT-FILM 
PROBE AND KULITE PRESSURE TRANSDUCERS 


Figures 18 and 19 show velocity traces created from 5.2429 seconds duration hot- 
film data samples at 90 percent speed. The data shown in Figures 18 was taken at a mass 
flow rate between peak efficiency and stall at 6.5 kg/sec. Figure 19 shows data that were 
taken near stall at a mass flow rate of 5.75 kg/sec. 
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Figure 18. Instantaneous velocity trace of hot-film output at 90 percent speed between peak 

efficiency and stall 
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Figure 19. Velocity trace of hot-film data at 90 percent speed near stall 


A power spectrum was created using a fast Fourier transform (FFT) of the hot- 
film data mentioned previously and the results can be seen in Figures 20 and 21. The 
power spectrum for the case between peak efficiency and stall shows a drop off in power 
level with increasing frequency, typical of free stream turbulence measurements. At a 
frequency of 8,938 Hz a spike occurred which corresponded to the blade-passing 
frequency (22 blades at 24,376 RPM). 
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Figure 20. Power speetrum of hot-film data at 90 pereent speed between stall and peak 

efficieney 

In the FFT of the near-stall condition the same overall shape appeared with the 
spike at the blade-passing frequency. Flowever there was slight evidence of the once-per- 
rev frequency of 406 FIz as well as an event at half this frequency, as annotated on Figure 
21 . 
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Figure 21. Power speetrum of hot-film data at 90 percent speed near stall 

An autocorrelation was conducted on the data taken at 90 percent speed at a mass 
flow rate between peak efficiency and stall. A plot of the autocorrelation is shown in 
Figure 22, which shows an overall exponential drop off of the signal with time. This too 
is typical of free stream turbulence. The high frequency oscillation superimposed on the 
exponential decay is once again the blade-passing frequency. The area under the curve is 
an indication of the size of the length scale of the inlet turbulence. 
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Figure 22. Autocorrelation of hot-film data at 90 percent speed between stall and peak 

efficiency 

The autocorrelation conducted on the hot-film anemometry data taken near stall at 
90 percent speed is shown in Figure 22. Besides the same overall trend as in Figure 21, a 
low frequency oscillation is apparent. This oscillation corresponded to a frequency of 
200Hz which was about 50 percent of the rotor speed. This was convincing evidence of 
the presence of a modal stall precursor oscillation. Previous measurements by Gannon et 
al. (Ref 5), showed that the subsequent stall cell rotated at 60 percent of the rotor speed. 
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Figure 23. Autocorrelation of hot-film data at 90 percent speed near stall 

The precursor oscillation was also found in the Kulite data near stall. A fast 
Fourier transform of Kulite data was used to create a contour plot of the power spectrum 
over time. Figure 24 shows data taken at peak efficiency. While the once-per-rev and 
blade-passing frequencies (or log thereof) were evident, the precursor frequency was not 
evident. 
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Figure 24. Power spectrum contour plot of Kulite data at 90 percent speed at peak efficiency 


Figure 25 represents data taken near stall at 90 percent speed. The precursor 
oscillation, once per revolution, and the blade-passing frequency, can be seen in this 
contour plot. 
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Figure 25. Power speetrum eontour plot of Kulite data at 90 pereent speed near stall 


Figure 26 is a eontour plot of Kulite data taken at 95 pereent speed and elose to 
stall. These data were taken at a mass flow rate that was not as elose to the stall margin 
as the data represented in Figure 25, but preeursor oscillation was also observed. Flot- 
film data at 95 percent speed are given in Appendix D. The measurement duration for 
these data samples was too short and yielded inconclusive frequency domain results. 
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Figure 26. Power spectrum contour plot of Kulite data at 95 percent speed near stall 


D. STEAM-INDUCED STALL AT 70% SPEED 

Initial testing of steam-induced stall was conducted at 90 percent speed with the 
steam pipe connected directly to the inlet duct that connected the settling chamber to the 
test compressor. This method yielded unsatisfactory results due to the fact that the steam 
was bypassing the throttle and increasing the mass flow rate through the compressor. 
The additional mass flow prevented stall rather than inducing stall in the compressor 
rotor. This method was abandoned when it was discovered that stall was easily achieved 
by venting steam outside and allowing it to be ingested through the throttle. The mass 
flow rate of the steam ingested during this initial attempt was not measured but the data 
from the steam-induced stall are presented in Appendix C. 

Figure 27 shows the hot-fdm, raw-voltage signal recorded by the high speed data 
acquisition system from the IFA 100, going into steam-induced stall at 70 percent speed. 
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Figure 27. Raw-voltage hot-film signal going into steam-indueed stall at 70 pereent speed 

Figure 28 is the raw-voltage of the Kulite signal going into steam-indueed stall at 
70 pereent speed. Beeause both the Kulite and hot-film data were reeeived by the VXI 
mainframe, the times seales on Figures 27 and 28 are identieal. 
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Figure 28. Raw-voltage Kulite signal going into steam indueed stall at 70 percent speed 
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A fast Fourier transform was used to create a power spectrum of the hot-film data 
as seen in Figure 29. Blade-passing frequency, once per revolution frequency and the 
rotating stall cell frequency were all identifiable from this trace. The corresponding raw 
signal of these data is given in Appendix E. 



Figure 29. Power spectrum of hot-fdm data during steam-induced stall at 70 percent speed 


The fast Fourier transform of the hot-film data was then evaluated over time to 
create a waterfall plot as seen in Figure 30. Figure 31 is the equivalent waterfall plot 
created from the Kulite data on the same time scale as the hot-film waterfall plot. 
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Figure 30. Waterfall plot of hot-film data going into steam-induced stall at 70 percent speed 
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Figure 31. Waterfall plot of Kulite data going into steam induced stall at 70 percent speed 
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Contour plot of hot-film data going into steam-induced stall at 70 percent speed 
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Figure 33. Contour plot of Kulite data going into steam-induced stall at 70 percent speed 
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Contour plots were also ereated from the fast Fourier transform. Figure 32 and 33 
show the hot-film and Kulite data respeetively. 

The ehange in speed of the eompressor going into steam-indueed stall is shown in 
Figure 34. 



Figure 34. Change in eompressor speed during steam-indueed stall 


The transient pressure measured in the steam line was used to ealculate the mass 
flow rate of the ingested steam. Figure 35 shows the pressure ehange in the steam pipe as 
well as the temperature ehange in the inlet of the eompressor over time during the steam- 
indueed stall experiment. A pressure of 475 kPa was reaehed in the steam pipe prior to 
releasing steam into the eompressor. By evaluating the ehange in pressure with respeet to 
time during the first seeond, a steam mass flow rate of 0.06 kg/see into the eompressor 
was ealeulated. 
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Figure 35. Steam pressure and inlet temperature ehange during steam-induced stall at 70 

percent speed 


A steam-induced stall of the rotor was observed at 70 percent speed at a mass 
flow rate of 4.30 kg/sec. Stall was not observed at higher mass flow rates using the given 
amount of steam. Transient data taken during steam ingestion but at a higher mass flow- 
rate which did not induce stall are given in Appendix F. Figure 36 shows the compressor 
map with the single point at 70 percent speed near stall that was measured prior to the 
steam ingestion experiment. At that point the compressor rotor stalled, hence the 

reduction in stall margin defined as ^peay_^\ calculated to be 

from 0.30 to 0.25, i.e. a 17 percent reduction in stall margin 
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Figure 36. Compressor map with 70 percent speed steam-induced stall margin 
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VI. CONCLUSIONS 


Hot-film measurements were suecessfully eonducted for the first time in the inlet 
to the transonic compressor rig over the rotor speed range of 70 to 95 percent of design 
speed. An average turbulence intensity of 2.26 percent with a variance of 0.287 was 
determined in the inlet flow of the compressor. Although variations in the turbulence 
intensity were observed, these variations were small and did not show a correlation with 
mass flow rate nor speed. 

Performance measurements were carried out at 94 percent speed from open 
throttle to near stall. These data were taken to establish a speed line closer to 100 percent 
speed typical of an engine at cruise condition. These data follow those previously taken 
with respect to the trend of the total pressure ratio, and the value of the measured peak 
efficiency. 

A precursor to stall oscillation was observed in the hot-film anemometry data as 
well as Kulite pressure transducer data. The precursor oscillation appeared to be at about 
half the rotor speed, and was only observable when the compressor was near stall. 

The rotor was successfully stalled at 70 percent speed as a result of steam 
ingestion. The release of steam was controlled and the mass flow was determined. 
Synchronized transient data, from the hot-film probe and Kulite pressure transducer, were 
recorded prior to and during the steam-induced stall event. Stall was not observed at 
higher mass flow rates using the given amount of steam and therefore the effect of the 
steam on the stall margin was determined. 
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APPENDIX A: HOT-FILM ANEMOMETRY CALIBRATION 

PROCEDURE 


IFA 100 Setup 

1. Record current temperature (Ta) and pressure (Pa). 

2. Insert shorting probe and measure cable resistance (Re). 

a. Pres [RES MEAS] 

b. Zero with [OPERATE RES] control knob 

c. Pres [RES MEAS] 

d. Record Rc then press [ENTER] 

3. Replace shorting probe with hot wire probe 

4. Repeat step #2 but DO NOT PRESS ENTER 

5. Press [OPERATE RES] and adjust to operating resistance on box. 

6. Press [BRIDGE COMP] and set with [BRIDGE] control knob according to the 
following table: 


Sensor Type 

Standard 1 

TE5 

35 

-P2 

83 

-PI1.5 

115 

-PI5 

120 

Metal Clad 

775 

-lOA 

70 

-10 

70 

-20 

115 

-60 

250 


7. Press [RUN] 

8. Start wind tunnel 

9. Turn [CABLE COMP] knob counterclockwise until OSC light comes on then 
turn clockwise until it goes off plus another half turn. 

note: it is recommended that step 9 be performed with the wind tunnel 
velocity closer to the upper bound of expected experimental 
velocity. 

10. The Low Flow Voltage (Eo) and High Flow Voltage (Em) need to be acquired in 
order to calculate GAIN and OFFSET. Eo and Em are the lower and upper 
bounds of the velocities expected during the experiment. Run the tunnel at both 
of these velocities and record the voltages. 

11. On the IFA main screen select [CALIBRATION] then select [PROBE DATA] 

a. Select [OPEN CAL FILE]. This creates a calibration file. It is 
recommended that this file be named after serial number of probe. 

b. Select [GAIN & OFFSET] 

c. Enter the voltages recorded in step 10 in the high flow and low flow box. 

d. Select [CALCULATE] and record gain and offset values. 
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12. In the CAL METHOD block, select option 2, “Acquire E and Type dP.” 

13. Select [CAEIBRATE] 

14. Enter current conditions on the left hand side of the screen remembering to select 
an appropriate dP unit. The dP is the difference between static and stagnation 
pressure and is used by the software to calculate the average velocity for each 
calibration point. 

15. On the right hand side of the screen, enter the number of calibration points that 
will be taken. 

16. Eor each calibration point, enter dP and reading and then press enter or click 
mouse outside of dP box. The lEA software should have calculated a velocity the 
dP reading. Press [Acquire] to record calibration point. 

17. After all calibration points have been taken, close graph and press [Next Screen]. 

18. Select [Curves]. 

19. At bottom of screen, in the “fit” block, select King’s law. Record King’s law 
coefficients in calibration sheet for future data analysis. 
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Hot Wire Probe Calibration Data Sheet 




Kings Law Coefficients 
A= -26.45895 


B=23.24134 


n=0.10000 
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Hot Wire Probe Calibration Data Sheet 



Ta 

Ambient 

Temp 

Pa 

Atm Pressure 

Rc 

Cable 

Resistance 

Ro 

Res. w/ probe 

Rop 

Oper. Res 

Bridge 

Eo 

L Flow 

Em 

H Flow 
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APPENDIX B: MATLAB M-FILES 


The following file was used to ealculate the average velocity and turbulence 
intensity as well as plot the Fast Fourier Transform of the data. 

function [turbint,avevel]=hotwirefft(king,file) 

% This function will calculate the average velocity, 

% turbulance intensity, and the Fast Fourier Transform. 

% format for file is ' 03NOV5 . E0005 ' 

% kings law coeff must be in a matrix with probes as rows, (a b n) 

% for example k=[.87726 .31691 .52] 
close all 

a=dlmread(file); 
b=dlmread(file); 

%This program is writen for sample size of 4096 

for n=l:4096 

a(n,3) = ( ( (a(n,2) *a(n,2) ) -king (1,1))/ (king (1,2) ) ) '^ (1 / king (1,3) ) ; 

end 

avevel_l=mean(a(:, 3) ) 
sd_l=std(a(:, 3) ) ; 
turbint_l=(sd_l*100)/avevel_l 

Figure(1) 

plot(a(:,1),a (:,3) ) 
xlabel('sec') 
ylabel('m/sec') 

%the following performs the Fast Fourier Transform 
voltage=b(:,2); 

N=length(voltage) -1; 
deltat=.00001; 

Y = fft(voltage) ; 

Y(l) = [] ; 
n=length (Y); 

power = abs (Y (1 :n/2) ) .'^2; 
nyquist = .5; 
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freq = (1:N/2)/(N*deltat); 


Figure(2) 

loglog(freq,power) 

xlabel{'Hz') 

ylabel { ' (m/sec) '^2/Hz ' ) 

The following file was used to plot the Fast Fourier Transform and Auto 
Correlation of the data. 

auto_fft.m 

close all 
clear all 

a=dlmread{'23NOV05.E0006'); 
time=a (:,1); 
voltage=a {:, 2 ); 

%This will plot the raw voltage trace 
Figure(1) 

plot(time,voltage) 
xlabel('time') 
ylabel('voltage') 
title('Raw Signal') 

N=length(voltage)-1; 
deltat=.00001; 

Y = fft(voltage); 

Y{1) = [] ; 
n=length{Y); 

power = abs (Y (1 :n/2) ) .'^2; 
nyquist = .5; 

freq = {1:N/2)/{N*deltat); 

Figure (2) 

loglog(freq,power) 
xlabel('freq') 
ylabel('power') 
title('Periodogram') 

% this will plot the auto correlation of the data 
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Figure(3) 

auto_time=deltat*lags; 
plot(auto_time, C) 
grid on 


Find Loc 

% Find_loc 

% m-function file to find the position of the trigger for each 

function [Loc] = Find_loc(time,tach); 

% Number of samples 
samples = length(time); 

% Trigger level is calculated 
Trig = mean([min(tach) max(tach)]); 

% Location of trigger points and correction to exact trigger 
% timing point 

Loc = find( tach(2:samples)<Trig & tach(1:samples-1)>Trig ); 

% Location of time of start of rev 


Movave.m 

% m-function file to calculate the moving averages at a 
% particular point 
% x,y data 

% points ahead and behind central one ie 0 return same data, 1 = 
% 3 points, 2=5 points 

% n, polynomial to fit, 0 = average, 1 = linear, 2 = parabolic 
% etc 

function [y_avg] = mov_avg(x,y,points,N) 

y_avg = y; 

x_poly = zeros(1,(2*points+l)); 
y_poly = zeros(1, (2*points + l)); 

for i = (points+1):((length(x)-points)-1) 
for j = -pointsipoints 


once per rev 
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x_poly(j+points+1) = x(i+j); 
y_poly(j+points+1) = y(i+j); 

end 

y_avg(i) = mean([max(y_poly) min(y_poly)]); % Most effective method, just take 
the mean of the max and min 

end 

% Leading points 
if points > 0 

% Leading few points 
for j = 1: {2*points + l) 
y_poly(j) = y (j); 

end 

% Leading moving average 
for i = lipoints 

y_avg(i) = mean([max(y_poly) min(y_poly)]); 

end 

% for i = lipoints 

% Trailing few points 
for j = {2*points+l):-l:l 

y_poly{j) = y{length(x)-j+1); 

end 

% Trailing moving average 

for i = (length(x)-points):length(x) 

y_avg(i) = mean([max(y_poly) min(y_poly)]); 
end % for i = lipoints 
end % if points > 0 

plotdata.m 

% m-file to plot a certain time part of the stall data file 

function [fred] = Plot_data(fig_no,Raw_data) 

% Time period is defined 

time_start = Raw_data(1,1); % Start time of sample 

time_end = Raw_data(end,1); % End time of sample 
%time_start =16.3 % User defined start time 

%time_end =17 % User defined end time 


% Time period entry points are found 
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temp = find(Raw_data(:,1)>time_start & Raw_data(:,1)<time_end); 

% Figure is plotted and trimmed 
%if bool_new_fig 

Figure(fig_no); close; Figure(fig_no); 

%else 

%Figure(fig_no) ; 

%end 

%if subplot_no ~= 0 

% subplot(subplot_tot ,1, subplot_no) 

%end % if subplot_no ~= 0 

plot(Raw_data(temp,1),Raw_data(temp,3),'b') 

%hold on 

xlabel('Time [s]'); ylabel('Raw Voltage signal [V]') 
grid on 

temp = axis; 

%h = line ( [time(Loc(1 : end) , Data_column) 

time(Loc(1 : end) , Data_column)] ', (ones(size(Loc(1 : end)) )*( [temp(3) temp(4 ) ] ))'); 

%set (h,' Color ',[ 0 00]); % Makes colour of line black 

%h = line([time(Loc(2 : end) , Data_column)-time_err time(Loc(2 : end) , Data_column)- 
time_err] ', (ones(size(Loc(2 : end))) *( [temp(3) temp(4)])) ' ) ; 

%set (h,' Color ', [0 00]); % Makes colour of line black 

%axis([temp(1 ) temp(2) 0 0.7]) 

%if subplot_no ~= 0 
% axis([0 120 0.5 1.3] ) 

%end % if subplot_no ~= 0 


plotdataHW.m 

% m-file to plot a certain time part of the stall data file 


function [fred] = Plot_data(l 

% Time period is defined 
time_start = Raw_data(1,1); 
time_end = Raw_data(end,1); 
%time_start = 16.3 
%time_end = 17 


g_no,Raw_data) 

% Start time of sample 
% End time of sample 
% User defined start time 
% User defined end time 
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% Time period entry points are found 

temp = find(Raw_data(:,1)>time_start & Raw_data(:,1)<time_end); 

% Figure is plotted and trimmed 
%if bool_new_fig 

Figure(fig_no); close; Figure(fig_no); 

%else 

%Figure(fig_no) ; 

%end 

%if subplot_no ~= 0 

% subplot ( subplot_tot ,1, subplot_no) 

%end % if subplot_no ~= 0 

plot(Raw_data(temp,1),Raw_data(temp,6),'b') 

%hold on 

xlabel('Time [s]'); ylabel('Raw Voltage signal [V]') 
grid on 

temp = axis; 

%h = line ( [time(Loo(1 : end) , Data_column) 

time(L og ( 1 : end) , Data_column)] ', (ones(size(Loo(1 : end))) *( [temp(3) temp(4)] )) ') ; 

%set (h, ' Color ', [0 0 0]); % Makes colour of line black 

%h = line ( [time(Loc(2 : end) , Data_column)-time_err time(Loc(2 : end) , Data_column)- 
time_err] ', (ones(size(Loc(2 : end))) * ([temp(3 ) temp(4)])) ' ) ; 

%set (h, ' Color ',[ 0 0 0]); % Makes colour of line black 

%axis ( [temp ( 1) temp (2) 0 0.7]) 

%if subplot_no ~= 0 
% axis ( [0 120 0.5 1.3]) 

%end % if subplot_no ~= 0 

Plotrpm.m 

% mfile to plot RPM through the stall 

function [fred] = Plot_rpm(fig_no,RPM) 

%temp = find(Raw_data (:, 9)>0) ; 

Figure(fig_no);% close; Figure(fig_no); 

%plot(Raw_data (temp, 1) , Raw_data (temp, 9)) 

%plot(RPM( :, 1 ), RPM( :, 2 ) , ' k ') ; 
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plot(RPM( \, 2 ) , ' k ' ) ; 
hold on 

xlabel('No. of Revolutions'); ylabel('Fraction of 90% Speed'); title('RPM through 
stall at 90% speed') 


Processdata.m 

% m function file to organise the data of the Raw data file 

function [time,volts,tach,Loc,RPM] = Process_data(Raw_data) 

% Data is sorted into simpler to use groups 

time = Raw_data(:,1); % Kulite time data 

volts = Raw_data(:,2:5); % Kulite raw voltage data 

tach = Raw_data(:,7); % Kulite trigger voltage 

samples = length(time); % Number of samples 

% Trigger level is calculated 
Trig = mean([min(tach) max(tach)]); 

% Location of trigger points and correction to exact trigger 
% timing point 

Loc = find( tach(2 : samples)<Trig & tach(1 : samples-1)>Trig ); %location of time of 
start of rev 

Hz = (length(Loc)-1)/(time(Loc(end))-time(Loc(1))); 

% Frequency of rotor revolution over the sample period 

% If Loc (location) if at the beggining of the sample it is 
% discarded 
if Loc(l) < 3 

Loc = Loc(2:end); 

end 

% Last trigger is discarded to ensure trailing zeros resulting 
% from probe lining up do not effect the calculations. 

Loc = Loc(l:(end-1)) ; 

% Timing correction to ensure that each sample begins at the 
% correct time 

d_time_Loc( :, 1) = time(Loc+1)-time(Loc) ; 

% Time interval between trigger and next time interval 
d_tach_Loc (:, 1) = tach(Loc+1)-tach(Loc) ; 

% Ramp slope between trigger and next time interval 
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m 


= d_tach_Loc./d_time_Loc; % Slope 


c 


= tach (Log); 


% Intercept 


time_err 


(Trig - c)./m; 


% Error in trigger 


% timing 


% Error in trigger timing is converted to be directly subtracted 
% from the period 

% time_err = (time_err(2 : end)-time_err(1 : end-1 ) ) ; 

% The RPM based on each trigger is calculated 
RPM = time(Loc(2:end)) - time(Loc{1:end-1)); 

RPM = RPM + (time_err(2:end)-time_err(1:end-1)); 

RPM = 60./RPM; 

%Hz_old = Hz; 

RPM = [time(Loc(2 : end)) RPM]; 

% First column is time signal and the second is the RPM 

% Data is moving averaged a few times to remove the wiggle 
%for i = 1 : N_mov_avg 

% [Hz(:;.2)] = mov_avg (Hz ( : ;. 1) , Hz ( : ;. 2 ) , 11) ; 

%end % for i = 1 :, N_mov_avg 

% [Hz ( : , 2) ] = mov_avg (Hz(:;.l),Hz(:;.2);.l;.l); 

% [Hz ( : , 2 ) ] = mov_avg (Hz ( : , 1) , Hz ( : , 2) ;. 11 ) ; 

% [Hz ( : , 2 ) ] = mov_avg (Hz ( : ;. 1) , Hz ( : , 2 ) ;. 11) ; 

% This is to correct the RPM 
%time_err = time_err + (Hz( :, 2)-Hz_old) ; 

% [ 1 . /Hz_old 1./Hz(:;.2) 1./ (time (Loc (2 : end) ) -time (Loc (1 : end-1) ) +time_err) ] 

% Data is converted to Hz or RPM 
%Hz {:, 2 ) = 1 . /Hz (:. 2) ; 


Stall 70.m 


%Stall90 

%M-file to pull in stall data and plot it out 
% 70% speed at surge 

close all; 
clear all; 
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%Raw data 


tic 

Raw_data=dlmread{'Dx2005_l123_0 958_7 0_6_10.csv' 
toe 

%tach=Raw_data( :, 8) ; 

%time=Raw_data( :, 1) ; 
data=Raw_data(:,2:6); 


'A6..H400000'); 


% Constants 
RPM_design = 27085; 

RPM_Hz = RPM_design/60; 

[time,volts,tach, Log, RPM] = 
[RPM{:,2) ] 


% Design RPM in RPM 
% Design RPM in Hz 

Process_data(Raw_data); 
mov_avg{RPM(: , 1) , RPM(: , 2) , 1, 1) ; 


%Input value for number of revs 

B=input('Enter a binary minus one (such as 511 for one rev) to select the number 
of revs:'); 


timeshortl=time(1:B, 1) ; 

deltatl=(timeshortl(end,1)-timeshortl(1,1))/(length(timeshortl)); 
%deltatl=(timeshortl(end,1)-timeshortl (1, 1))/(length(timeshortl)-1) ; 
[Loc] = Find_loc(time,tach); 
kulitel=data{:,1); 


%Figure 1; 

%_A particular column of the data is plotted 
plot_data(1,Raw_data); 
title{'Stage Surge 90%'); 


%Figure 2; 

%_The RPM is plotted 
plot_rpm(2,RPM/RPM_design) ; 
title{'Stage Surge 90%'); 


%Figure 3; 

%The beginning of stall is graphed on a waterfall plot 
%Figure(3) ; 
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count=l; 

Al=200; 

A2=600; 
for A=A1:A2; 

%_Locates triggers 
LocA= {Log (A): {Log (A) +B)); 

yfftA=kulitel {LogA, 1) ; 

Ya=fft(yfftA); 

Ya(l) = [] ; 

powerl=abs{Ya{1:{length(yfftA) / 2 ) ) ) .^ 2 } 
Nl=length(yfftA); 
nl={l:(Nl) 12 ) 

freql={nl)/{(Nl)*deltatl); 

timel=(time (Log (A)))*(ones(size(freql))) 

if 0 

hold on; 

Figure(A); 

semi logy(freql,power1, 'r' ) ; 


end 

freq_fall(Gount,:)=freql'; 
power_fall(GOunt,:)=powerl'; 
time_fall(GOunt,:)=timel'; 

GOunt=GOunt+l; 


end 

power_fall=loglO(power_fall); 
freq_fall =loglO{freq_fall); 

%_Graphs a waterfall plot 

%waterfall(freq_fall,time_fall,power_fall) ; 

%_Graph settings 
ylabel('log(Frequeney)'); 
xlabel('Time'); 
zlabel('log(Power)'); 
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%title('stage Surge 70%: Entering Stall at Triggers 10-30'); 
grid on; 


%Figure 4 

%The beginning of stall is graphed on a contour plot 

%Figure (4); 

count=l; 

%for A=A1:A2; 

% contourf(time_fall,freq_fall,power_fall,15); 
shading flat 

%_Graph settings 

xlabel('Time') ; 

ylabel('log(Frequency) ') ; 

%title{'Stage Surge 70%: Entering Stall at Triggers 10-30'); 
grid on; 

%axis{[time(Loc(Al)) time(Loc{A2)) 2 5]); 

count=count+l; 

%end 
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APPENDIX C: 90 PERCENT AND SPEED STEAM-INDUCED 
STALL: HOT-FILM AND PRESSURE DATA 

Hot-film file 10NOV05.E0015 


Raw Signal 




freq 
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Autocorrelation Autocorrelation 



Tlme(s) 



Time(s) 
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Waterfall plot of 10NOV05.E0015 hot-film data file 


4.5 



2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5 2.55 2.6 


Time 


Waterfall plot of Kulite data. 



1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 

Time (sec) 
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vdtage 


APPENDIX D: 95 PERCENT SPEED NEAR STALL: HOT-FILM 

DATA 


Data file 03NOV05.E0010 


Raw Signal 



a 


10 


10 


10 


10 ' 


10 


10 


10 


10 


10 

freq 



10 


10 
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Autocorrelation Autocorrelation 



Tlme(s) 



Tlme(s) 
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APPENDIX E: 70 PERCENT SPEED STEAM INDUCED 
STALL: HOT-FILM DATA 


Hot-film velocity trace from the IFA 100 and processed using MATLAB program 
from Appendix B 
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power voltage 


APPENDIX F: 70 PERCENT SPEED STEAM INGESTION 
WITHOUT STALL: HOT-FILM DATA 


Raw Signal 



0 1 2 3 4 5 6 


time 



freq 
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Autocorrelation Autocorrelation 



-- 

-- 

-- 

-- 

-- 

— 






















_ 1 







_ 1 

_ 

__ , 

_ , 

_ 1 

_ 


0.8 


0.6 


0.4 


0.2 


0.5 


1 1.5 

Tlme(s) 


2.5 


1 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 



- 0.01 - 0.005 0 0.005 0.01 0.015 0.02 


Tlme(s) 
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